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Abstract 

Following Refs. [1, 2j, in the Littlest Higgs model (LHM), the heavy photon is supposed to 
^_ be a possible dark matter candidate and its relic abundance of the heavy photon is estimated 

^^ , in terms of the Boltzman-Lee- Weinberg time-evolution equation. The effects of the T-parity 

^D \ violation is also considered. Our calculations show that when Higgs mass Mh taken to be 300 

^^ ' GeV and don't consider T-parity violation, only two narrow ranges 133 < Majj < 135 GeV 

and 167 < Mai^ < 169 GeV are tolerable with the current astrophysical observation and if 
135 < Majj < 167 GeV, there must at least exist another species of heavy particle contributing 
to the cold dark matter. As long as the T-parity can be violated, the heavy photon can decay 
into regular standard model particles and would affect the dark matter abundance in the uni- 
verse, we discuss the constraint on the T-parity violation parameter based on the present data. 
_H _ Direct detection prospects are also discussed in some detail. 

Qh' Keywords Little Higgs dark matter, time evolution. 

^ ! I. Introduction 

^ . The Littlest Higgs model |3J has been proposed for solving the hierarchy problem of Stan- 

'sj" ! dard Model (SM) recently. Unfortunately, the original Littlest Higgs model still suffers from 

r^ ' severe constraints from precision electroweak fits and the fine-tuning of the Higgs boson mass 

^T) • is necessary. To avoid the difficulty, a discrete Z2 symmetry named as "T-parity" |4j (just as 

R-parity in SUSY) is introduced to the model. This symmetry in the little Higgs models en- 
dow SM-like particles T-even parity while all heavy gauge bosons and scalar triplets are T-odd. 
trr* ■ Therefore, it is completely natural to see that the problem in the model can be avoided. An 

interesting feature of the Littlest Higgs Model with T-parity is that the T-odd particles need to 
be pair-produced and will cascade down to the lightest T-odd particle (LTP) which means that 
^ ' the LTP is guaranteed to be stable. At the same time, one of the most fundamental problems 

in cosmology and particle physics today is what is the nature of the Dark Matter. A flood of 
discussions point out that the Dark Matter should have the following characters - non-luminous, 
non-baryonic, non-relativistic, and electrically neutral [5| [6l 171 [HI [9| [TO] . 

The heavy photon in the Littlest Higgs model indeed meets all criteria, thus it can be an 
ideal candidate of cold dark matter (CDM). The relic abundance of the dark matter in the 
Littlest Higgs model with T-parity in the thermal relic scenario has already been evaluated by 
Hubisz and Birkedal [HIS]. To be consistent with the observation on CDM, the mass of of heavy 
photon is estimated around a few hundred GeV. 

In this work, following Refs.fU [2], we re-evaluate the abundance of heavy photons in our 
universe by means of the Boltzman-Lee- Weinberg equation where the T-parity violation effects 
are taken into account. As a reasonable consideration we set the mass of Ah in the range: 100 
GeV < Majj < 300 GeV . We find that a rather narrow window for the mass of heavy photon 
meet the requirements of the observation, outside of which other kinds of dark matter particles 
are necessary. The paper is organized as follows: after this introduction, we formulate the time 
evolution of the dark matter in the Littlest Higgs model with and without T-parity violation in 
section H and then we present our numerical results. Direct detection prospects are discussed 
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in section III. The section IV is devoted to our brief conclusion and discussion. 



II. The time evolution of the dark matter 



The time evolution of the dark matter is described by the Boltzman-Lee- Weinberg time- 
evolution equation [11], which is the basis of study on the abundance of dark matter for a 
possible candidate particle (heavy photon in this work) at the present stage: 
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where < av > is the average value of the heavy photon pair-annihilation cross section times the 
relative velocity;R is the cosmic scale factor and R/ R is the Hubble constant H; n is the number 
density of cold dark matter in the thermal bath and no is the equilibrium density: 
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the expression of no is different to our former paper [12], where neutralino is supposed to be the 
cold dark matter particle which is a fermion. Instead, Ah we employ here is a boson and obeys 
the Boson-Einstein statistics. 

Obviously, in eq.([T|), the T-parity is conserved, so that the heavy photon as the lightest 
T-odd particle does not decay. The evolution equation should be modified if we consider the 
T-parity violation effects: 
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where T = 1/r is the decay width of the Ajj which decays into SM particles via T-parity viola- 
tion. The < av > means the thermally averaged annihilation cross section. 

In this work we only consider that the dark matter particles annihilate into the SM particles 
in s-wave, thus under the non-relativistic limit, the cross section av is simply given as 

av = av\^^^, +av\^^ +av\HH + <7v\tt . (4) 



Due to relatively large uncertainties of the present data, we restrict ourselves to only eval- 
uate the cross sections of the 2< — > 2 processes which would be overwhelming over the 2 to 3 
modes, and the corresponding Feynman diagrams are shown in Fig. [TJ 




Figure 1: Feynman diagrams for the annihilation of the dark matter (Ah) in the 2 <-> 2 pro- 
cesses. Diagrams (1) to (4) are those which give the largest contributions to the annihilation 
< av > for the ranges of Mah and Mh that we examine. Diagram (5) contribution to the 
annihilation coefficient < av > is much smaller than the other four, so we ignore it in numerical 
computation [2]. 
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Table 1: The interaction vertices Littlest Higgs model that appear in this paper. 



For the readers' convenience, we present the vertices which we are using in Table , where 
the coupling g = 2-^/7ra/cos^^. We can directly obtain av of the diagrams shown in Fig. [H 
The subscript "CM" here denotes the center-of-mass frame. 
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where Hi = Mf /M\ and Mah is the mass of Ah, 9^ is the SM Weinberg angle, a keeps the 
meaning that it is in SM and AIh is the mass of SM Higgs boson and Th is the width of the 
Higgs. If Mah > ^TT't, the annihilation of heavy photon pair into a pair of top quarks is possible, 
and 
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If Mah > Mh, it also can annihilate into pairs of Higgs bosons, with the following cross 
section [2j, 
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where the subscript "CM" denotes the center-of-mass. The dependence of av on the masses 
of Higgs and heavy photon is depicted in Fig. [21 
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Figure 2: This plot depicts the dependence of the < av > on the mass of Higgs and heavy 
photon. From the lightest to the darkest regions, the magnitude oi < av > decreases as shown 
in Fig. 



In the numerical computations, we take a few typical values for mass of Higgs as M/f =240 
GeV and Mh=300 GeV, as shown in Fig. O , 
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Figure 3: the magnitude of < ctv > with Mj:/=240 GeV (left) and Mh=300 GeV (right). 

It is noted that we have not considered the contributions of the T-parity violation effects to 
< (TV > yet. Later in this work, we will show that to be consistent with the observational data 
on dark matter density in our universe, if the heavy photon is the only constituent of the dark 
matter, the T-parity violation cannot be large. Thus when we evaluate < av >, the diagrams 
related to the T-parity violation can be safely neglected. There is a sharp peak in the figure 
which is due to the resonance effect with 4M^ — M^ ~ 0. 

For readers' convenience, let us briefly repeat the procedure for solving Eq. (3). Following 
the literature, one needs to rewrite the first term by means of the equation of state, according 
to the Friedmann-Lemaitre equations 
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where H(t) is the Hubble parameter, p is the energy density and A is the cosmological constant. 
Energy conservation Tf^'^ = leads to another useful equation. 



p = -3H{p + p) . 
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This equation can be written as p = —3(1 + w)pR/R and is easily integrated to yield 

p ex i?-3(i+«') , (11) 

where the state parameter w = p/p is 1/3 at the early universe dominated by radiation and 
when matter dominates. 

For the radiation-dominant universe, one has 
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where g* = g^ + -gp is the effective number of degrees of freedom. The g^ and g^ denote the 
degrees of freedom for boson and fermion respectively. 

After the universe converted from radiation-dominance to matter-dominance, the relation 
RT = constant still holds [11], and 



p = f5T^ 
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where constant [3 can be fixed by the present energy density 

Pr = f3T^ = 0.mPc , 

where T^ = 2.725 K [IB] is the present microwave background radiation temperature, ^mh'^ 
0.127 is the present-day density parameter for pressureless matter and the critical density pc 

——- = 1.05 X IQ-^h^ GeV-cm-3 [13], thus 
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We have the modified evolution equation for the radiation-dominant stage as 
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and the other equation for the matter-dominant universe 
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In Fig. U we plot the dependence of p/T'^ on temperature (the time). In order to illustrate 
our numerical results, without losing generality we set the mass of Higgs as Mh = 300 GeV and 
will make more discussions in the last section. 
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Figure 4: p/T^ vs T. where two states of radiation and matter are compared. Here Mh = 300 
GeV and Ma„ = 170 GeV. In case I, as an approximation, we use radiation-dominance equation 
(15). The matter-dominance equation (16) is used in case II when temperature blow 10^ K. 



The early universe was hot, dense and dominated by radiation, but as the temperature 
gradually decreased, non-relativistic matter eventually dominates the universe. Thus Eq.(15) 
and Eq.(16) describe different stages of the temperature evolution and the solutions of the two 
equations should be smoothly connected at a medium temperature of about lO'^ K. With this 
restriction we can obtain our numerical results which are shown in Figs. |3]~ Figs. [71 

By Fig. m one can notice that when temperature exceeds 10^^ K, f{x) (or p/T^) approaches 
to original fo{x) (or po))iiamely the curve almost does not change until the temperature drops 
to the decoupling temperature Tj. As the temperature decreases continually the curve is flat 
again for a long range until the temperature Ty where T-parity violation applies, then the curve 
drops abruptly. 

In this figure, for understanding the physics, we plot two curves. The line I is a solution of 
Eq. (15) where one assumes the universe is dominated by the radiation all the time, whereas 
the line II corresponds to the case where the two solutions of Eq.(15) and Eq.(16) are smoothly 
connected, namely the two stages of time evolution are properly dealt with. One can notice 
that without taking into account of the T-parity violation, line I and line II are almost overlap, 
but when one considers the T-parity violation, the line I drops earlier than the line II. It may 
influence the time evolution of dark matter in the future (billion year later, we suppose). 

The Feynman diagrams which are responsible for the heavy photon decaying into SM par- 
ticles are shown in Fig. [U where the vertex is T-parity violated. 




Figure 5: The Feynman diagrams for the decay of the heavy photon Ah- Diagrams (2) to (3) 
can be ignored since Zq, W± and H will decay furthermore and the phase space is too narrow. 
/ in diagram (1) stands for all quarks, leptons except the top quark as long as Ah is set as 170 
GeV in the numerical computations. 



The decay amplitude of diagram (1) of Fig. [5] is 
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where K, Pi, P2 denote the momenta of Ah, f and / respectively and Xf is the T-parity 
violation parameter. The decay width of Ah that decays to a pair of // is 
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Because ruf <C M^^ , we can ignore it safely in the expression and obtain 
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where Cf = 3 for quarks and cj = 1 for leptons. Finally, we get the total decay width of Ah 



Ttotai = Y.^f = 0.0796Ma^(5A2 + 2X1) = 0.0796Ma„X^ 



(20) 



where Ai, A2 denote the T-parity violation parameters of quarks and leptons respectively. In 
general, Ai and A2 may be different, but the equation shows that only their combination applies 
in practice as an effective parameter A. 

Because of T-parity violation, the curves dip down at certain low temperature that depend 
on the mass of the heavy photon and the T-parity violation parameter A. For an illustration, 
we choose three typical values of lifetimes of heavy photon which should be sufficiently long 
compared to the age of our universe, otherwise it may contradict to the observation of the dark 
matter density, as lO"^^, 10^^ ^nd 10"^° GeV"^ (about 2 x 10^^ 2 x 10^° and 2 x 10*^ years) cor- 
responding to the dimensionless T-parity violation parameter A = 3.00 x 10~^^, 3.00 x 10~^^, 
3.00 X 10-21 for Ma„ = 140 GeV; A = 2.72 x 10-^3^ 2.72 x 10-^2^ 2.72 x lO^^i fo^. ^^^ ^ 179 
GeV; A = 2.08 x lO'^^, 2.08 x 10-^2^ 2.08 x lO^^i fo^. j^^^ ^ 290 GeV respectively. The three 
cases are shown in Fig. [6l 
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Figure 6: p/T^ vs T. Different values for the heavy photon mass Majj and T-parity violation 
parameter A are taken. 



By measuring the anisotropy, the density of CDM is determined as |13j 
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where h is the Hubble constant in units of 100 km/(s-Mpc) and we can get the range of present 
CDM density pcDM, 
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49.30 eV • cm"-' • K^'^ < pcDhi/Tr < 58.12 eV • cm"^ • K 
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Obviously if the heavy photon is the unique constituent of CDM, it must account for the 
whole CDM energy density, and the curves in Fig. [U] must pass through or intersect with the 
range Eq. (21). 

The figure shows that there is only a narrow window for its mass to be allowed by the 
measurements on the dark matter density. For clarity, in Fig. \7\ we show the relation between 
the density of CDM, (p/T^) in the flat region of curves in Fig. [Hand the mass of heavy photon 
M^if The space between the two horizontal lines corresponds to the observational data with 
error tolerance. It is shown that without T-parity violation, i.e. heavy photon is absolutely 
stable, only two regions 133 < Ma^ < 135 GeV and 167 < Ma^ < 169 GeV can be consistent 



with astrophysical data, if 135 < Ma^ < 167 GeV, there must be at least another kind of heavy 
particles contributing to the rest CDM density. It is interesting to notice that as Ma^ changes, 
the cross section < va > varies and near the resonance shown in Fig. 3, it is proportional to 
(M^ /r|) and becomes very large. Back to the master equation Eq.(3), the solution of po/T^ 
can reach almost zero near the resonance. Therefore the two allowed ranges of M^^ are sepa- 
rated and the dependence of po/T'^ on Ma^ is not monotonic. 
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Figure 7: Relation between the CDM density p/T^ = n x Majj/T'^ and the heavy photon mass 
Ma^ ■ The middle region is allowed by the present astronomical observation. 

If T-parity violation indeed exists, the situation would be different. For 1 a error tolerance, 
nearly all values of Ma^ except 135 < Ma^j < 167 GeV is possible, because Ma^j can decay via 
T-parity violation into light SM particles. However, on the other aspect, its relic density must 
pass through the allowed region which is ranged by the solid lines (1 a) or dashed lines (3 a) 
of Fig. [7] , that is constrained by the present astronomical observation data. This as well sets a 
restriction to the T-parity violation parameter A for different Ma^- 




Figure 8: The relation of Ma„ and F as Ma„ >170 GeV. The shaded area is allowed by 
the present observation, with here we only consider the situation of Mai^ > 170 GeV for an 
illustration in the 3a error tolerance. 



To be explicit, we present the relation between Ma^j and the decay width A which is con- 
strained by data, in Fig. 8. Unfortunately, the relation cannot be written in a simple analytic 
form, so that we illustrate it graphically. As Majj gets larger, the lifetime should be shorter, 
i.e. the T-parity violation effect should be stronger (the violation parameter is larger) accord- 
ingly. Thus if we assume Ma„ to be a component of CDM, the T-parity violation parameter is 
rigorously constrained as it decays slowly to meet the present observation. 

It is interesting to investigate how the initial conditions influence the evolution, in other 
words, at very high temperature, if we change the initial density of heavy photon, whether we 
can obtain different relic density at present time. We choose three different initial values at high 
temperatures and draw the evolution curves in Fig. [9l the numerical results indicate that the 
evolution is almost independent of the initial conditions. It implies that the relic density is only 
related to the decoupling temperature and mass of the CDM particles. 




Figure 9: p/T^ on T with various initial conditions. Here Ma^ is taken as 170 GeV. 



III. Direct detection experiment 



There are many attempts to detect dark matter, either directly or indirectly. In this work, 
we only concern the possibility of directly detecting the dark matter flux. 

In direct detection experiments one measures the dark matter flux via elastic scatterings 
of the heavy photon off nuclei (proton) . In this work, we would theoretically evaluate the cross 
section in the Littlest Higgs model and estimate probability of observing such processes. 

We suppose that the heavy photon is the only component of the CDM flux, and then whether 
we can have sufficient data to find a clear signal depends on the flux density and the dispersive 
velocity of Ah in the solar neighborhood and the coupling of heavy photon with proton. This 
information about the concerned parameters is obtained by studying the relic density of CDM 
and that procedure was carried out in previous sections, we then evaluate the rate of events 
expected in an earth experiment (i.e. LTP-proton scattering events) per unit time, per unit 
detector material mass. The events number is approximately given by 

An = iV < cr > nUAte, (23) 



where N is the number of target protons in the detector, n = pq/Majj { Po = 0.3 GeV/crn? 
) [13] is the local heavy photon number density and < av >~< a >I;is the cross section for 
the scattering which is averaged over the relative velocity of Ah in the lab frame, v =< v^ >2 
is the dispersive velocity of dark-matter particles, At is the time of the detection, e is the de- 
tection efficiency. In order to get a plausible numerical result, we select 1000 kg water, thus 
A'" = NAW^g/{18g/mol) x 10, the local energy density of dark matter Ah{po) is 0.3 GeV/cm^ 
[13] . the mass of Ah{Majj), is 175 GeV and At is set to be one year. 
The formula of the average cross section < a > is shown as below: 



< a{v,) >= / f{vi)a{vi)d''vu (24) 



where vi = \vi\ is the velocity of the Ajj in the lab frame where the target proton is at rest, and 
14] is the local Maxwellian velocity distribution, by v =< v^ >2= 270 km/s [14^ [15] , we know 



v^ = 5.40789 X 10-^c2. 

The expression for the cross section is a bit more complicated 

Mh 4pf - M^ 

a{vi) = 5r(-ui)[arctan(---^) - arctan( S.)] ^ (25) 

i_ff 1/f 

where we take the Higgs mass Mh = 300 GeV, T^ is the decay width of Higgs, pi = \pi\ = 
MpVi/{MAH + ^p) is the momentum of Ah in the CM frame and 

9{vi) = / " , . (26) 



192^pIMhTh{^pI + Ml^ + ^pI + M2)2 
The Feynman diagram of heavy photon scattering off proton is shown in Fig. [TOl 



H 



Figure 10: Feynman diagrams for the directly detecting the dark matter {Ah)- 

Through equations (23)- (26), we know that if the detection efficiency e is about 0.1, during 
one year period, in 1000 kg water, there can only be 5.8 x 10~^^ signals. If we want to obtain 
58 signals, we need at least 10^^ tons water, that is to say, we need to build a 10 m high, 10^ m 
wide and 10^ m long pool, which sounds unreasonable. At the same time. Unfortunately, due 
to the extremely small interaction rates and the small energy transferred to the recoiled nucleus 
that if the dark matter particle is heavy photon, it is almost impossible to distinguish the events 
from backgrounds. 
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IV. Conclusion 

The purpose of measuring dark matter flux is not only to measure the flux itself, but also to 
explore new physics beyond the standard model. Indeed, as indicated in all literature, the CDM 
particles cannot be SM ones, therefore study on the CDM flux may provide valuable information 
about the new physics. It is lucky that the LHC will begin running in 2008, one may expect 
to find signals about new physics beyond the SM from a great amount of available data. Thus 
the cosmology constraint on the new physics will stand as complimentary to the accelerator 
experiments. 

There have been many proposals about the CDM candidates in literature, such as the SUSY 
particles neutralino [IS], sneutrino[I71[l8l[l9], axino[20l[2ll[22], technicolor-particles [Ml [Ml EH], 
heavy photon in the Littlest Higgs model [26l [271 [M] aiid even unparticles[29l[3Q], it is hoped to 
find an efficient way to distinguish signals originating from different models. In fact, by solving 
the Boltzman-Lee- Weinberg equation and making the results to be consistent with the obser- 
vational data, one can set constraints on the model parameters, such as the mass of the CDM 
candidate particle and its coupling to the SM particles if its decay is allowed. This estimation 
will provide information for the accelerator experiments which are searching for such particles 
on earth. Moreover, the probability of directly detecting of the dark matter particles on earth is 
another way to distinguish the various models. Namely, even though all the the models can be 
consistent with the astronomical observation on the relic dark matter, they may have completely 
different probabilities to be caught by the detector on earth. By measuring the dark matter 
fiux, we can rule out some possible candidates and then combining the data of LHC, one may 
eventually determine the identity of the new physics particle for CDM. 

In this work, a possible Dark Matter candidate - heavy photon in the Littlest Higgs model 
with T-parity has been discussed. We observe that different Higgs and heavy photon masses 
can result in totally distinctive average value of the heavy photon pair-annihilation cross sec- 
tion, which determines the time evolution of our universe. By employing a sample Higgs mass 
Mh = 300 GeV without losing generality, we present the relationship between the density of 
cold dark matter p/T^ and the temperature T for different stages. Thanks to the newly avail- 
able data and observation of astrophysics, we can set a stronger constraint on the mass of heavy 
photon. Only two ranges 133 < Majj < 135 GeV and 167 < Majj < 169 GeV are allowed 
without T-parity violation , whereas there must be at least another kind of heavy particle con- 
tributing to the rest CDM energy density for the range 135 < Majj < 167 GeV. Outside the 
three ranges, the Littlest model cannot meet the data. If dark matter particles are found to 
fall in the regions beyond 132 to 170 GeV (of course with some errors), then the Littlest Higgs 
model without T-parity violation would be ruled out by the cosmology data. If T-parity vio- 
lation was considered, the data would set rigorous constraints on T-parity violation parameter 
A. We also calculate the possibility of directly detecting of this type of Dark Matter, i.e. heavy 
photon on the earth. Our conclusion is that it is impossible to observe any signal due to small 
interaction rates and severe backgrounds. In other words, if heavy photon is the CDM particle, 
it can easily evade our detection on earth. Since different CDM particles have different proba- 
bilities to be detected by earth detectors, our results offer a way to distinguish the Littlest Higgs 
model from others. If our detectors on earth do not catch any dark matter fiux particles, the Lit- 
tlest Higgs model may survive, otherwise, one can erase such a model from our model collections. 
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